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ABSTRACT 

Recent results by the WMAP satellite suggest that the intergalactic medium was 
significantly reionized at redshifts as high as z ~ 17. At this early epoch, the first 
ionizing sources likely appeared in the shallow potential wells of mini-halos with virial 
temperatures Tvir < lO^K. Once such an ionizing source turns off, its surrounding HII 
region Compton cools and recombines. Nonetheless, we show that the "fossil" HII 
regions left behind remain at high adiabats, prohibiting gas accretion and cooling 
in subsequent generations of mini-halos. Thus, early star formation is self-limiting. 
We quantify this effect to show that star formation in mini-halos cannot account for 
the bulk of the electron scattering opacity measured by WMAP, which must be due 
to more massive objects. We argue that gas entropy, rather than IGM metallicity, 
regulates the evolution of the global ionizing cmissivity, and impedes full reionization 
until lower redshifts. We discuss several important consequences of this early entropy 
floor for reionization. It reduces gas clumping, curtailing the required photon budget 
for reionization. An entropy floor also prevents H2 formation and cooling, due to 
reduced gas densities: it greatly enhances feedback from UV photodissociation of H2. 
An early X-ray background would also furnish an entropy floor to the entire IGM; 
thus, X-rays impede rather than enhance H2 formation. Future 21cm observations 
may probe the topology of fossil HII regions. 



1 INTRODUCTION 



Tlie high optical deptli r = 0.17 ± 0.04 detected by the 
Wilkinson Microwave Anisotropy Probe ( WMAP) satellite 
has lent greater credence to the notion of an early period of 
star formation and reionization, Zr — 17 ± 8 (Kogut et al 
2003; Spergel et al 2003). If indeed the first stars formed 
at high redshift z ~ 20, they are expected to form in mini- 
halos^ with shallow potential wells, in which H2 coohng is 
dominant (Abel, Bryan & Norman 2000; Bromm, Coppi 
& Larson 2002). More massive halos with Tvir > 10*K, 
in which coUisional ionization and line cooling can operate, 
are expected to be very rare at these redshifts. Modulo the 
effects of UV feedback on H2 formation and coohng, stars 
forming in such halos could therefore play a dominant role in 
an early reionization epoch. A great deal of effort has gone 
into assessing the impact of UV feedback on H2 cooling, as 
well as the counter-vailing effects of positive feedback effects 
such as an early X-ray background (e.g., Haiman, Abel & 
Rees (2000); Ciardi, Ferrara & Abel (2000); Ricotti, Gnedin 
& ShuU (2002); Machacek, Bryan & Abel (2001)). 



^ For the purposes of this paper, a mini-halo is defined as any 
halo with T^jr < 10^ K which cannot cool via atomic line cooling. 



The main point of this paper is that the population 
of mini-halos is likely to be considerably sparser than pre- 
viously assumed. This is because mini-halo formation is 
strongly suppressed even inside the fossil HII regions of dead 
ionizing sources. Although such HII regions recombine and 
cool by Compton scattering with cosmic microwave back- 
ground (CMB) photons, they cannot cool to back to the tem- 
perature of the undisturbed intergalactic medium (IGM). 
Strong Jeans mass filtering takes place (Gnedin 2000), and 
subsequent mini-halos will no longer be able to accrete gas 
due to the smoothing effects of finite gas pressure. Thus, 
once any patch of the universe is ionized, it can no longer 
host any more mini-halos, even if it subsequently cools and 
recombines. In effect, the birth of the first stars leads to the 
demise of the mini-halo population: only one generation of 
stars can form within these shallow potential wells. 

While many authors have noted and commented on 
the Jeans mass filtering after full reionization (Bullock, 
Kravtsov & Weinberg 2000; Benson et al 2002; Somerville 
2002), the Jeans filtering in fossil HII regions after ionizing 
sources have turned off has not been studied. We argue that 
it will strongly suppress the mini-halo population, with the 
following interesting consequences: 
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• Impact of Mini— Halo Population on Reioniza- 
tion. Since the recombination time is shorter than the Hub- 
ble time at high redshift, tree ^ tn, and the ionizing 
sources, expected to be massive stars, have short lifetimes, 
this ~ 3 X lO^yr ^ tu, many generations of star formation 
are required to keep a given patch of the IGM ionized. How- 
ever, once the first generation of stars born in mini-halos 
dies out, subsequent generations will not be able to form 
in mini-halos in a previously reionized patch of IGM, even 
after the patch cools and recombines. This effect inevitably 
produces a non-monotonic reionization history with an early 
peak of partial reionization, followed by recombination and 
eventual full reionization. This is qualitatively similar to 
the reionization history derived by Cen (2003) and Wyithe 
& Loeb (2003), although the physical reason for the non- 
monotonic evolution in the present work is different (entropy 
injection, rather than a pop HI to pop II transition caused 
by a universal metallicity increase). Regardless of the extent 
of UV feedback effects on H2 production and cooling, the 
majority of stars which reionized the universe were hosted 
by more massive halos able to survive Jeans mass filtering, 
T^ir > few X lO^K. 

• Photon Budget for Reionization. Dense gas that 
collects in mini-halos can form a considerable sink of ioniz- 
ing photons by boosting the overall effective gas clumping 
factor Cii = {nl)/{ne)'^ (Haiman, Abel & Madau 2001; 
Barkana & Loeb 2002; Shapiro et al 2003). The clump- 
ing factor Cii increases rapidly at low redshift as structure 
formation progresses; thus, despite the higher mean gas den- 
sity n cx (1 + z)^ at high redshift, the recombination time 
tree = l/laUeCii) docs uot evolve strongly from z = 10 — 20. 
However, if the universe is filled up with fossil HII regions at 
high redshift, the photon budgets required to subsequently 
ionize it and to keep it ionized are much lower, since gas 
clumping in mini-halos is strongly suppressed. 

• Source Clustering. Since preheating boosts the 
threshold mass required for efficient gas cooling and star 
formation, it increases the mean bias of the early proto- 
galaxy population by a factor of a few (e.g., see figure 2 of 
Oh, Cooray & Kamionkowski (2003)). This is likely to in- 
crease the clustering amplitude of background fluctuations 
due to these faint unresolved early proto-galaxies, such as: 
the free-free background due to ionized halos (Oh 1999; Oh 
& Mack 2003), Sunyaev-Zeldovich fluctuations due to high 
redshift HII regions and supernovae winds (Knox et al 1999; 
Oh, Cooray & Kamionkowski 2003; Santos 2003), and the 
IR background due to stellar emission (Santos, Bromm & 
Kamionkowski 2002; Magliocchetti, Salvaterra & Ferrara 
2003; Haiman, Spergel & Turner 2003). The increase in 
clustering bias boosts the amplitude of rms fluctuations by 
a factor of several. This may be necessary, for instance, if 
thermal SZ fluctuations from high-redshift supernovae are 
to account for the small-scale CMB anisotropies observed 
by BIMA and CBI (Oh, Cooray & Kamionkowski 2003). 

• Star Formation History/Metallicity/IMF. Ini- 
tially, each mini-halo is expected to harbor a single mas- 
sive, metal-free star (no fragmentation is seen in numerical 
simulations by Abel et al. 2000 and Bromm et al. 2002) . Be- 
cause of the preheating, no further mini-halos can form in 



the entire comoving volume reionized by this mini-halo. As 
a result, the first generation of stars is expected to form as 
a population of single, isolated stars. Star formation ensues 
again only once more massive halos with deeper potential 
wells aggregate. These rare, high density peaks are likely 
to coincide with the highly biased regions where the first 
isolated stars had already formed. At the time of the for- 
mation of deeper potential wells, such sites are inevitably 
already polluted with metals, unless the first stars collapsed 
directly to black holes without associated metal production. 
As a result, the transition from pop III to pop II (metal free 
to normal) stellar populations is associated with halo mass 
scale. This is in contrast to scenarios (Cen 2003; Wyithe & 
Loeb 2003) where star formation can continually proceed in 
lower density peaks, which are far from the initial sites of 
star formation and still contain relatively pristine gas (al- 
lowing metal-free star formation to continue to relatively 
low redshifts). Given the factor of ~ 10 — 20 difference in 
the ionizing photon production efficiency per unit mass be- 
tween Pop HI and Pop II stars (Tumhnson & ShuU 2000; 
Bromm, Kudritzki & Loeb 2001; Schaerer 2002), this wiU 
have important consequences for the redshift evolution and 
topology of reionization. 

• 21cm Observations. It has been suggested that mini- 
halos will be observable at the redshifted 21cm line fre- 
quency both in emission (Iliev et al 2002a; Iliev et al 2002b) , 
and absorption (Furlanetto & Loeb 2002). If mini-halo for- 
mation is strongly suppressed, these two windows on small- 
scale structure during the cosmological Dark Ages will dis- 
appear. 21cm observations of mini-halos therefore provide 
an interesting probe of the topology of fossil HII regions: if 
mini-halos are seen, that comoving patch of the IGM has 
never been ionized. This is complementary to other observa- 
tions such as Gunn-Peterson absorption, which only probe 
the instantaneous ionization state. 

The rest of this paper is organized as follows. In 
we compute the "entropy fioor" due to early mini-halos. In 
S|3 we calculate the effects of this entropy floor on the gas 
density profile in mini-halos. In 21 discuss and quan- 
tify several important effects of this entropy fioor for reion- 
ization, and in S|3 we explicitly calculate the effect on the 
reionization history. We summarize our findings and discuss 
their implications in |H] Throughout this paper, we adopt 
the background cosmological parameters as measured by the 
WMAP experiment (Spergel et al. 2003, Tables 1 and 2), 
Q.m = 0.29, Q.A = 0.71, fit = 0.047, h = 0.72 and an ini- 
tial matter power spectrum P{k) oc fc" with n = 0.99 and 
normalization erg = 0.9. 



2 ENTROPY FLOOR DUE TO EARLY 
REIONIZATION 

2.1 Entropy floor in Fossil HII regions 

Early reionization introduces an 'entropy fioor' in the inter- 
galactic medium which impedes gas accretion and cooling in 
halos with Tvir < 10*K which cannot excite atomic cooling. 
Early reionization even introduces Jeans smoothing effects 
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Figure 1. The parameter dependence of the final temperature of 
a parcel of isochorically cooling gas. The solid lines indicate the 
results from the full non-equilibrium chemistry code, while the 
points indicate the approximate analytic solution from equations 
1331 and 1351 which is in excellent agreement. Our fiducial patch 
has an initial temperature Ti = 20, OOOK, cools at z = 19, for 
a time interval t = O.Gtuiz = 19) = 1.9 X lO^yr (corresponding 
to roughly Sz ^ 7), and lies at an overdensity 5 = 1. The gas 
cools most effectively at high redshift (when the Compton cooling 
time is short) and at low overdensity (when it decouples from the 
CMB at late times). For z=19 and S = 10, (which is perhaps 
characteristic of gas being accreted onto halos) the gas remains 
at ~ few X lO^K. 



in halos with Tvir > 10'*K, since the strength of the accretion 
shock is weaker in preheated gas, and a deeper potential well 
is required to coUisionally ionize the gas. The overall effect 
is to introduce a substantial core in the gas density profile. 
The effect is similar to the presumed preheating of the IGM 
by galactic winds or AGN outflows at low redshift. There, 
the finite entropy of the gas introduces a core in group and 
cluster gas density profiles and is thought to be responsible 
for the deviation from self-similarity in the observed cluster 
Lx — Tx relation (for a recent review, see Rosati, Borgani 
& Norman (2002)). Entropy is more fundamental than gas 
temperature in low-redshift clusters and high-redshift mini- 
halos: in both cases, entropy is conserved, since icooi ^ tn- 
We exploit this analogy, and deliberately employ language 
and techniques from the well-studied low redshift case to 
study high-redshift mini-halos. Here, we calculate the ex- 
pected level of the 'entropy floor' due to early reionization. 

Reionization is likely to be a highly stochastic process 
where sources ionize a patch of the IGM and then fade; the 
fossil HII region subsequently cools and recombines. What 
is the final temperature a fossil HII region can cool to? At 
the high redshifts of interest, Compton cooling off the CMB 
is by far the dominant source of gas cooling in the low den- 
sity IGM (e.g., at 2 = 19 and T = 10*K, tc ~ O.ltcooiS, 



where 5 is the gas overdensity, and tcooi is the atomic cool- 
ing timescale). The final gas temperature is therefore de- 
termined by the competition between Compton cooling and 
hydrogen recombination: after the gas recombines, it 'de- 
couples' from the CMB, and Compton cooling is no longer 
efficient. The Compton cooling timescale is independent of 
density and temperature: 

tc = 3mec(4a:eaTar^MB)"' = 1.4x 10^ {^^) ^ ""-^^^^ 
whereas the recombination time 

f,ec = 3.9 X 10' (-^) ' 5~^3^e\r (2) 

is shorter for higher density gas. Thus, denser gas 'decou- 
ples' from the CMB earlier and freezes out at higher tem- 
peratures. For the redshifts of interest, Ucc(5 = 1) > tc, 
allowing for substantial cooling below ~ lO'^K. 

We can study the temperature evolution of a cooling 
and recombining parcel of gas T{t,Ti, 5, z) as a function of 
time t (or redshift interval Sz), initial temperature Ti, gas 
overdensity Si and redshift z. In Figure 2 we show the de- 
pendence of the final temperature on these parameters by 
solving the full set of chemical evolution equations for pri- 
mordial gas. We neglect H2 chemistry as a trace UV flux is 
sufficient to photodissociate H2 at these low IGM densities. 
Our fiducial patch has an initial temperature Ti = 20, OOOK, 
cools at z — 19, for a time interval t = 0.6ti^(z = 19) (cor- 
responding to roughly Sz ~ 7), and lies at an overdensity 
5 = 1; each panel shows the dependence on one of the vari- 
ables holding the others fixed. We do not take into account 
time- dependent variation of the gas density arising from in- 
fall onto non-linear structures or Hubble expansion, which 
would provide adiabatic heating/cooling of the gas; however 
the effects of this may be roughly estimated by consider- 
ing some mean overdensity 5 of the gas during its evolution. 
The final temperature is almost independent of T, as long 
as t ^ tc and the gas was initially fully ionized. However, 
the final temperature declines at higher redshifts and low 
overdensities S, since ic/iroc oc (1 + z)^^ S , and the gas cools 
faster than it recombines. 

In Appendix I, we develop an analytic expression for 
T{t, Ti, 5, z) taking into account only Compton cooling and 
recombination which is remarkably accurate; it is given by 
the points in Figure [H This allows one to estimate the final 
temperature quickly without evolving the coupled differen- 
tial equations. It deviates from the true answer significantly 
only at high overdensity and/or low redshift, when atomic 
line cooling becomes comparable to Compton cooling. In any 
case, for metal-free gas, line cooling alone can only cool the 
gas down to ~ 5000K. The only processes which could alter 
the final entropy of the gas significantly are H2 cooling and 
metal line cooling. These are more efficient than Compton 
cooling only in regions of high overdensity, in the halo core, 
and so do not affect the initial entropy floor. 

Once cooling becomes ineffective, the gas evolves adi- 
abatically. In accordance with convention, we shall refer to 
the quantity: 
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as the 'entropy' of the gas, even though it is not strictly 
the thermodynamic entropy 5* oc log(K). This is a useful 
convention because K is conserved when the gas evolves 
adiabatically. Thus, it is conserved during Hubble expan- 
sion, as well as during accretion onto halos (provided the 
accretion shock is weak, gas will be accreted isentropically) . 
In Figure |5| we show the dependence of the final entropy 
K on redshift (the dependence on Ti,t follow simply from 
Figure 0. The gas has significantly lower entropy at higher 
redshift, since Compton cooling is more efficient. Ifowever, 
the final entropy at a given redshift depends only weakly on 
the overdensity S. Denser gas remains at higher tempera- 
ture, since it recombines faster than it cools, and the com- 
bination of higher temperature and higher density roughly 
cancel K oc T^-^^^ We can therefore ignore the weak 5 de- 
pendence of K, and assume that it depends only on z. The 
entropy floor is thus roughly independent of the details of 
structure formation and gas density distribution. 

It is useful to define a parameter which compares 
the IGM entropy floor to the entropy generated by grav- 
itational shock heating alone. Let us define the quan- 
tity k = Kigm/Ko, where Ko = Tvir/n^rvii)^^^ , and 
n(rvir) = i^b/i^7n)pNFw{rvii)/itJ,mp) (and pnfw is the 
NFW (Navarro, Frenk & White 1997) dark matter den- 
sity profile). This is the entropy due to shock heating alone 
at the virial radius; the justification will become clearer in 
the ensuing section. As K increases, the Jeans smoothing 
effects due to finite IGM entropy become increasingly more 
significant. This parameter K{Kigm ,Tvir, z) will be used 
extensively in the following sections, and it is useful here to 
get a sense of what values of K are expected. In Figure |21 
we plot K for a Tvir ~ 9000K halo as a function of redshift, 
using Kigm{z) shown in Figure |21 Since K oc T^^^, and a 
Tvir ~ 9000K halo is about the most massive that would still 
not experience atomic cooling, the solid line depicts a lower 
limit on K for all mini-halos. Appropriate values for K for 
smaller halos can be read off simply from the K oc T~^^ scal- 
ing. We see that for most halos, K > 1; although Kigm{z) 
falls at high redshift, this is somewhat offset by the fact that 
typical potential wells are much shallower, and thus Ko also 
falls rapidly (as seen by the dashed lines, which depict K 
for 2a and 3a fluctuations). 

Typical values for K are illustrated further in Figure 
13] where we show the mass weighted fraction of mini-halos 
which have entropy parameters less than a given K, given 
by: 

J^^f^dMMjdn/dM) 

/(< K, z) = —j- ; 4 

/^"" dMM{dn/dM) 

where M{K) is the mass corresponding to K for a given 
Kiom{z), Mj{z) is the cosmological Jeans mass (e.g., see 
Barkana & Loeb (2002)), Mu{z) is the mass corresponding 
to a Tvir = 10*K halo, and dn/dM is the Press-Schechter 
mass function. Virtually all halos have if > 1 in fossil HII 
regions at all redshifts of interest {z <, 20), and median 
values of K are much higher. We shall soon see that such 
halos are subject to very substantial Jeans smoothing effects. 
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Figure 2. Dependence of the final entropy K = T jr?!^ on red- 
shift z and overdensity h. Points correspond to the analytic solu- 
tion from equations I33l and lij5l As in Figure our fiducial patch 
has Tj = 1, <5 = 1, ^ = 19, t = 0.54/^(2 = 19). The dependence on 
Ti and t are the same as in Figure^ The gas retains more entropy 
at lower redshift, since Compton cooling is less efficient. However, 
the final entropy shows only a weak dependence on overdensity <5: 
the gas cools out at higher temperatures at higher densities, and 
the increased temperature and density roughly cancel. 



2.2 Entropy Floor due to X-rays 

Reionization by X-rays (Oh 2001; Venketesan et al 2001) 
would produce a warm (few x 100 — lOOOK), weakly ionized 
IGM with an entropy floor similar to that in fossil HII re- 
gions. Such X-rays could arise from supernovae, AGN, or 
X-ray binaries, or in more exotic models with decaying mas- 
sive sterile neutrinos (Hansen & Haiman 2003). The universe 
is optically thick to all photons with energies: 

E < Sthick = 1.8 (iy^) xKi'keY (5) 

where xm is the mean neutral fraction, and we have as- 
sumed cr^ oc Thus, all energy radiated below iSthick 
will be absorbed. The relative efficiency of UV photons and 
X-rays in setting an entropy floor deserves detailed sepa- 
rate study; here we discuss some salient points. UV photons 
are an 'energetically extravagant' means of producing an en- 
tropy floor. Most of the energy injected by UV photons is 
lost to recombinations and Compton cooling at high redshift; 
we see in Figure that the gas typically Compton cools to 
Tfloor ~ few X lOOK at the redshifts of interest. Thus, only 
~ 10~^ of the heat injected by UV photons is eventually 
utilized in setting the entropy floor; all energy expended in 
heating the gas above Tfloor is 'wasted' (it is dumped into 
the CMB, where it may eventually be observable as a spec- 
tral distortion, Fixsen & Mather 2002). On the other hand, 
in the weakly ionized gas produced by X-rays, both the re- 
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Figure 3. The evolution of the dimensionless entropy parameter 
K = KiQi^/Ko as a function of redshift, assuming the gas cools 
at the mean IGM density 5 = 1. Here Ko is the gas entropy at the 
virial radius due to shock heating alone, and ifiGM is the entropy 
of the fossil HII regions as computed in Figure 1^1 The solid line 
describes a Tvir = 9000K halo, the most massive in which atomic 
cooling is still not important. It therefore defines a lower limit on 
K. The entropy parameter can be simply scaled {K oc TT^ ) for 
halos with lower virial temperatures. The dashed lines show K for 
2(7 and Scr halos at that redshift; the vast majority of mini-halos 
have k > I. 



combination time and the Compton cooling time are longer 
than the Hubble time. In particular. 



tc = 1 



l + z 
15 



tH 



(6) 



and almost none of the entropy injected is lost to cooling. 
Furthermore, a larger fraction of the X-ray energy goes to- 
ward heating rather than ionization (in general, fewxO.l of 
the energy of the hot photo-electron created by an X-ray 
goes toward heating; this fraction quickly rises toward unity 
as the medium becomes progressively more ionized (ShuU 
& van Steenberg 1985)). Thus, if eboi(X - ray)/eboi(UV) > 
lO"'^ (where eboi is the comoving emissivity), X-rays could 
be comparable or even more effective than UV photons in 
setting an entropy floor. The relative emissivities of UV and 
X-rays is unknown, but for supernovae could be as high as 
(Oh 2001): 



<;boi(X - ray) 
eboi(UV) 



0.1 



(7) 



where /osc is the escape fraction of ionizing UV photons 
from the host halo, and fx is the fraction of supernova ex- 
plosion energy which goes into soft X-rays, either through 
inverse Compton scattering of CMB photons by relativistic 
electrons (Oh 2001), or free-free emission from the hot SN 
remnant. 



Figure 4. The mass weighted fraction of mini— halos which have 
entropy parameters < K, for z = 10,15,20, assuming JsTigmC^) 
shown in Figure]^ Virtually all halos have K > 1, and the median 
value of K is substantially higher. Thus, the vast majority of halos 
accrete gas isentropically. 



To make a quick estimate, let us (fairly conservatively) 
assume that fx ~ 3% of the explosion energy of a super- 
nova goes into soft X-rays. Of this, /heat ~ 50% of the energy 
goes into heating; the rest goes into secondary ionizations 
and atomic excitations. A supernova releases Ez ~ O.SMeV 
in explosion energy per metal baryon, relatively independent 
of metallicity (a Pop III 'hypernova' produces ~ 100 times 
more energy, but also ~ 100 times more metals than a stan- 
dard type II SN). X-ray heating thus results in an entropy 
floor KiGM ~ (/x£'z/hoat^)/T^^/^ or: 



KiGM ~ 20eVcm 



Vo.osy V 



/he 



0.5 



10-3Zg 



(8) 



where Z is the mean metallicity of the universe. This is 
comparable to the entropy of fossil HII regions, but has a 
much larger filling factor, of order unity. A metallicity of 
Z ~ 10~"^Zq corresponds to roughly ~ 1 ionizing photon 
per baryon in the universe, both for Pop II and Pop HI 
stars. Due to recombinations, the filling factor of ionized re- 
gions will be of course considerably less than unity. However, 
the filling factor of fossil HII regions, which is the relevant 
quantity, will also be less than unity by a factor /overlap, 
which represents the overlap of HII regions with fossil HII 
regions which have recombined. Energy spent in reionizing 
fossil HII regions is 'wasted' in terms of establishing an en- 
tropy fioor. This factor is likely to be large because early 
galaxy formation is highly biased, and higher mass halos 
(which can resist feedback effects) will be born in regions 
already pre-reionized by earlier generations of mini-halos. 
A fossil HII filling factor of order unity with little overlap 
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can only be achieved if the formation of mini-halos is highly 
synchronized (see discussion in section (SJ . 

The mean free path of X-ray photons generally exceeds 
the mean separation between sources, becoming comparable 
to the Hubble length for ~ 2keV photons. The entire uni- 
verse is thus exposed to a fairly uniform X-ray background, 
and the entire IGM acquires a uniform entropy floor, with an 
amplitude which scales with the amount of star formation, 
as in equation (|SJ. Even regions far from sites of star forma- 
tions, which have never been engulfed in an HIT region, will 
be affected. By contrast, in pure UV reionization scenarios, 
there are large spatial fluctuations in entropy, which depend 
on the topology of reionization, and the redshift at which a 
comoving patch was last ionized. 

Because of the relative uncertainty of the amplitude of 
the X-ray background, we use the entropy floor associated 
with fossil HII regions in the rest of this paper. The mass 
fraction of affected mini-halos scales with the flUing factor of 
fossil HII regions. This is therefore a minimal estimate; the 
filling factor could approach unity if X-rays are important. 

We now consider the effects of a finite entropy floor on 
mini-halo gas density profiles. 



3 GAS DENSITY PROFILES 

Once Compton cooling and radiative cooling become ineffi- 
cient, the gas evolves adiabatically. We can therefore com- 
pute static equilibrium density profiles, and see how they 
change as a function of the entropy floor. The models we 
construct are in the spirit of Voit et al (2002); Tozzi & 
Norman (2001); Oh & Benson (2003); Babul et al (2002), 
which match observations of low-redshift cluster X-ray pro- 
files well. 

Naively, one might assume that only gas which remains 
at temperatures comparable to the virial temperatures of 
mini-halos would suffer appreciable Jeans smoothing. Since 
in many cases the IGM can cool down to few x lOOK, one 
might assume that this level of preheating would have negli- 
gible effects on the density profile of gas in mini-halos. This 
is false: the important quantity is not the temperature but 
the entropy of the gas. Since gas in the IGM is heated at 
low density, it has comparatively high entropy. Gas at mean 
density which is heated to temperatures: 

^-->^o(3=)(25o)"''^ 

(where 5 is the overdensity of the gas in the mini-halo in the 
absence of preheating) will have entropy in excess of that ac- 
quired by gravitational shock heating alone. Its temperature 
will therefore exceed the virial temperature after infall and 
adiabatic compression. As the level of preheating increases, 
gas at progressively larger radii in the halo undergoes Jeans 
smoothing effects. We now calculate this in detail. 

We first construct the default entropy profile of the gas 
without preheating. We assume that in the absence of heat- 
ing or cooling processes, the gas distribution traces that of 
the dark matter, an ansatz which is indeed observed in nu- 
merical simulations (e.g., (Frenk et al 1999)) (the dark mat- 
ter is assumed to follow the NFW (Navarro, Frenk & White 



1997) profile). This assumption becomes inaccurate at the 
very center of the halo, where finite gas pressure causes the 
gas distribution to be more flattened and less cuspy than 
the dark matter density distribution. In particular, even in 
the absence of preheating the IGM has a flnite temperature 
after decoupling from the CMB and cooling adiabatically: 
Tmin(z) « 2.73(1 + Zd)[(l + z)/{l + Zd)f, where the matter- 
radiation decoupling redshift Zd ~ 150. This gives rise to a 
finite entropy floor: 

K^i,, = 4.6 X 10~^eVcm^ (10) 

independent of redshift. Thus, there will be a finite core in 
the gas density profile even in the absence of preheating. 
In regions where gas traces the dark matter, hydrostatic 
equilibrium gives the entropy profile due to shock heating 
as: 

Kshoak{r) = JJ-^ X (11) 

■ r _ GM(.KW ,,^P^,^r^,^.(,^^^)\ 

The final entropy proflle is therefore K{r) = 
max[A'min, -K'shock('')]. Note that the mean molecular 
weight IS n — 0.59 for fully ionized gas and ^ = 1.22 
for fully neutral primordial gas; in the paper we are 
dealing with the case where the ionization fraction is small 
and therefore use the latter figure. This results in virial 
temperatures which are higher by a factor ~ 2. 

What is the effect of preheating on the entropy pro- 
file? If the infall velocity does not exceed the local sound 
speed, then the gas is accreted adiabatically and no shock 
occurs; the gas entropy A'igm is therefore conserved. If gas 
infall is supersonic, then the gas is shocked to the entropy 
^shock computed in equation (IIH . Tozzi & Norman (2001) 
found that the transition between the adiabatic accretion 
and shock heating regime is very sharp. To a very good 
approximation K(M) = max(7('shock(Af), -ffiGM); with this 
new entropy profile we can compute density and tempera- 
ture profiles. In fact, for most high redshift mini-halos the 
'preheating' entropy exceeds the shock entropy even at the 
virial radius, so the gas mass is accreted isentropically. This 
occurs when the entropy fioor exceeds: 

-'(5^)(4)"'(^)"'.'- 

where the gas overdensity 5 ~ 50 at the virial radius in 
the absence of preheating is a weak function of the NFW 
concentration parameter c. Although c depends weakly on 
the collapse redshift, for simplicity we shall assume in this 
paper that c = 5 for all mini-halos. 

We have compared our entropy profiles calculated with 
equation (IIH with the prescription in Tozzi & Norman 
(2001) (and also adopted by Oh & Benson (2003)). In this 
method an accretion history for a halo is prescribed via ex- 
tended Press-Schechter theory; this allows one to compute 
the strength of the accretion shock and thus the gas en- 
tropy (using standard Rankine-Hugoniot jump conditions) 
for each Lagrangian mass shell. The two methods agree ex- 
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tremely well; we therefore use equation IIH for both speed 
and simplicity. 

Given an entropy profile, from hydrostatic equilibrium 
the density profile of the gas is then given by: 

p(r) = A:(r-)-'/^ X (13) 

p,.,,.,-ui^£-*-^*(/,-] 

where K = P/p^^^ = kB{fJ.mH)~^^^K. The temperature 
profile can then be determined from T = Kp1~^ . 

The solution of this equation requires a boundary 
condition which sets the overall normalization, here ex- 
pressed as P(rvir)- There is some ambiguity in this choice. 
The often used boundary condition A4g = fsMhaio at 
r = Tvir is unphysical as it does not take into account 
the suppression of accretion due to finite gas entropy. We 
make the following choice. Let us define Pshock(''vir) = 
(r2i,/nm)p(rvir)/(/i?7i£f )fcs7vir, the prcssure at the virial 
radius due to shock heating alone. For isTshockl^'vir) > 
^iGM, the final conditions at the virial radius are not 
strongly affected by the entropy floor, since the shock boosts 
the gas onto a new adiabat. We therefore set P(rvir) = 
fshock(?'vir). The bouudary condition must change when 
KiGM > -R'shock(''vir), whcu accretiou takes place isentrop- 
ically, and the entropy floor is fundamental in determining 
the gas pressure. In this case, P(rvir) (which is essentially 
a constant of integration) is chosen so that p(r) ^ p as 
r —> oo. The latter boundary condition implicitly assumes 
that hydrostatic equilibrium prevails beyond the virial ra- 
dius. This is questionable, but is arguably reasonable: the 
gas sound speed Cs = (l-ft'iGMP^"^)^''^ in preheated gas is 
much higher than in cold gas, which is essentially in free- 
fall. The sound-crossing lengthscale over which hydrostatic 
equilibrium can be established, Lsc ~ CstH, is: 




which is ~ 5 times larger than the virial radius for a ~ 6000K 
halo at the same redshift. Indeed, by deflnition Lsc > ^vir 
for Jeans smoothing effects to be important. 

Having established the appropriate boundary condi- 
tions, we can now compute detailed gas proflles. The en- 
tropy, density, pressure and temperature profiles as a func- 
tion of f = r/rvir are shown in Fig|S] These profiles are of 
course universal and independent of halo mass once K is set, 
if the weak dependence of the NFW concentration param- 
eter c with mass is ignored. As the entropy floor increases, 
the central pressure and density decline, while the central 
temperature increases. We see that reasonable values of the 
entropy floor K > 1 produce dramatic effects on the gas 
density profile of the halo, smoothing it out considerably. 

We now use these density profiles to compute how the 
accreted gas fraction fg = (Mg/Mhaio)/(n6/f^m) scales with 
the entropy parameter K. This is shown as the dark solid line 
in Fig|S| Again, because of the self-similarity of the problem, 
this plot is valid for mini-halos of all virial temperatures 
at all redshifts, provided K is appropriately re-scaled. It is 
reassuring to see that fg is continuous at iC" = 1, when we 
switch from one boundary condition to another. This need 
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Figure 5. The dimensionless entropy K = K/Ko, pressure 
P = P/Po, temperature T = T/To, and density p = p/po, as 
a function of radius f = r/ r^ir ■ Ko , Po,To are the values of 
these quantities at Tvir without preheating, while po = ph, the 
mean baryonic density. The entropy profile K(f) uniquely speci- 
fies P(f),T(f), p(f), independent of halo mass or redshift. 

not have been the case, and gives us confidence that we 
handle the transition to isentropic accretion correctly. We 
see that realistic levels of the entropy floor (as computed in 
f|5J causes a substantial depression in gas fractions in mini- 
halos. 

It is interesting to compare our derived gas frac- 
tions with other estimates. For the case where A'igm > 
^t^shock('"vir), accrction takes place isentropically. The halo 
will therefore accrete gas at roughly the adiabatic Bondi 
accretion rate (e.g., Balogh, Babul & Patton (1999)): 

Mt « IMtvXG^mLioK-^H (15) 

The total accreted gas mass is then roughly Adgas ~ 
mm{fMttH, {0.b/i^m)Mhaio), where / is some unknown nor- 
malization factor which takes into account the fact that the 
total halo mass is not constant but was lower in the past 
(and hence, that the gas accretion rate was lower in the 
past). Another estimate which is a good flt to the results of 
hydrodynamic simulations is (Gnedin 2000): 

'""^ [l + (2V3_l)Mi/2/M,„,„]3- ^"'> 

There is only one free parameter: Afi/2, the mass of the 
halo in which the gas mass fraction fg — 0.5. Gnedin (1999) 
shows that M1/2 is well approximated by the "filtering mass" 
Mf- However, Mf depends on the unknown thermal history 
of the IGM. To make a self-consistent comparison, we com- 
pute Ml/2 with the density profiles computed with our fidu- 
cial boundary conditions. Interestingly, we find that 
roughly corresponds to the halo mass when accretion begins 
to take place isentropically A" ~ 1. 



© 0000 RAS, MNRAS 000, 000-000 



8 Oh & Haiman 



Isentropic 
Accretion 



1000 - 





1000 ? 



Figure 6. The gas fraction within the virial radius /gas = 
(Mg/Mjiaio)/(f2i)/fim) as a function of the entropy parameter K. 
The solid line indicates our fiducial boundary conditions; the dot- 
ted line corresponds to the Gnedin (2000) fit to numerical simu- 
lations; while the dashed line corresponds to the Bondi accretion 
prediction (valid only for isentropic accretion, K > 1). All three 
show good agreement. The self-similarity of this problem implies 
that the computed gas fractions applies to mini-halos of all virial 
temperatures, provided K is appropriately scaled. For illustrative 
purposes, appropriate values of K for a mini-halo of Tvir = 5000K 
ai z = 20, 13 arc shown (see Figure 1^. 

The results are shown in Figure |S| AH three estimates 
agree well (note that the normahzation / of the Bondi ac- 
cretion prediction is a free parameter; the plot shown is for 
/ = 0.3). The widely used Gnedin (2000) fitting formula 
predicts even lower gas fractions (and as we will see, clump- 
ing factors) at high entropy levels K ^ 1. On the other 
hand, the slope of the fg (K) relation for our boundary condi- 
tions agrees very well with the Bondi accretion prediction in 
this regime. Our boundary conditions therefore yield fairly 
conservative estimates of the effects of preheating. More- 
over, unlike these other estimates, we are able to compute 
detailed density profiles, which is crucial for some of our 
later calculations. 

There are two other quantities which are of particular 
interest when computing density profiles. One is the central 
density pc, which affects the ability of mini-halos to form 
H2 in the face of UV photo-dissociation. Another is the gas 
clumping of the halo, defined as: 

Chalo ~ [ . \ ■ (17) 
W 

where the brackets indicate a volume averaged quantity, 

(X) = ^ f drATTr^X. (18) 
^ Jo 

Note that Chaio > 1 always. The clumping factor Chaio plays 



Figure 7. The effect of increasing the entropy parameter K on 
the central gas density pc = pel P (bottom panel), and gas clump- 
ing Cj^{fjj (top panel), as defined in equation 1171 . The solid lines 
correspond to our fiducial boundary conditions, while the dashed 
lines correspond to adjusting P(r^i,.) to reproduce the gas frac- 
tions from the Gnedin (2000) fit to numerical simulations (the 

results for the Bondi accretion prediction are almost identical to 
~ 1/2 

the solid curve). Both p^ and Cj^^^j^ decline rapidly as K increases. 



a central role in determining the photon budget required for 
reionization; we evaluate the global clumping factor in the 
following section. In Figure|7| we show the effect of increas- 
ing the entropy parameter K on the central density and 
clumping factor. Both decline rapidly with K. We shall use 
these two results in the following sections. 

We now use these gas density profiles to compute global 
effects of mini-halo suppression. 



4 GLOBAL EFFECTS OF MINI-HALO 
SUPPRESSION 

4.1 Suppression of Collapsed Gas Fraction and 
Gas Clumping 

An entropy fioor suppresses the fraction of gas which is 
bound within mini-halos. As we discuss in H4.4I the mean 
21cm emission from mini-halos is directly proportional to 
this global gas fraction; if it is strongly suppressed the sig- 
nal will be unobservable. The global collapsed gas fraction 
is: 

/gas = - / dM-^fth,i,,^{K)M (19) 

where /s — Q.i,/Q,m, fhaio,g{K) is the halo gas fraction as 
in Figure 1^1 and pt is the comoving baryon density. This is 
shown in the bottom panel of Figure El Curves are shown 
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Figure 8. The mass function of baryons, as in equation I2UI for 
z = 10, 15, 20. Dark solid curves show the mass function for 
^IGM = ^min, light Solid curves are for redshift-dependent en- 
tropy ^^10111(2) as in the top panel of Figure 2, and dashed curves 
are for a fixed entropy fioor iCiGM = 3eVcm^. Note that the 
baryonic mass function is very similar at all redshifts for Kiqm{z). 
For 

Tvir > lO^K, an entropy floor has no effect on the mass func- 
tion, which reverts back to the dark solid curves (we do not extend 
the light solid and dashed curves above the corresponding halo 
masses) . 



for no preheating Kigm = -ft'min, fixed values of the en- 
tropy floor ifiGM = 1, lOeVcm^ (which correspond to the 
IGM settling to some temperature at a given redshift and 
remaining at constant entropy thereafter), and the redshift 
dependent entropy T^igm (z) shown in the top panel of Fig. 
2, which reflects the increasing entropy of gas which Comp- 
ton cools at late epochs. For realistic values of Kigm, the 
collapsed gas fraction is suppressed by 1-2 orders of magni- 
tude. 

It is interesting to plot the effect of preheating on the 
baryonic mass function. This can be computed simply as: 

dN dN dMh ^20) 



dN_ dMh 

dMb ~ dMh ^ dMb ' 

We show this in Fig. |H1 at 2; = 10,15,20 for KiGmiz) as 
in the top panel of Fig. 2, and for a fixed entropy floor 
Kigm ~ 3eVcm^. As expected, the effect of preheating is 
most drastic at low masses. Note that the baryonic mass 
function is very similar at all redshifts for Kigm{z); this can 
be also seen in the bottom panel of El where the collapsed 
fraction in mini-halos does not change appreciably with red- 
shift. 

Another quantity of great interest is the global gas 
clumping factor C = (n^)/{n)^. Gas clumping shortens the 
recombination time tree ~ l/{anC) and thus increases the 
total number of photons per baryon required to achieve 
reionization. It is given by: 



where fv is the collapsed fraction by volume of mini-halos, 
CiGM « 1 is the clumping factor of the IGM, Chaio(^) 
is the halo clumping factor as calculated in the previ- 
ous section and shown in the top panel of Figure |7| and 
Vhaio = (47r/3)r^ir is the volume of a mini-halo. The factor 
(tevap/tff) deserves some explanation. Mini-halos only con- 
tribute to recombinations when they are photoionized. How- 
ever, once they are exposed to ionizing radiation, they will be 
evaporated, essentially in the sound crossing time of photo- 
ionized gas (Shapiro, Raga & Mellema 1998; Barkana & 
Loeb 1999). As in Haiman, Abel & Madau (2001), we there- 
fore set the evaporation time tovap ~ rvir/(10kms~^), and 
weight the halo clumping factor by the duty cycle tevap/in. 

The clumping factor in equation 1211 is a lower limit to 
the total gas clumping; there will be an additional contribu- 
tion from larger scales as well. However, the clumping due 
to minihalos is dominant at the high redshifts of interests. 

Our result is shown in the top panel of Figure|3] For the 
no entropy floor Kigm = Kmin case (top curve), gas clump- 
ing is much lower at high redshift, since fewer mini-halos 
have collapsed. Once a patch of IGM is reionized early at 
high redshift, then an entropy floor is established and sub- 
sequent gas clumping is suppressed. Thus, early star forma- 
tion reduces the total photon budget required to achieve full 
reionization. We see that gas clumping is reduced by an or- 
der of magnitude even for the low entropy levels ~ leV cm'^ 
associated with reionization at high redshift z ~ 20. For the 
redshift-dependent entropy Kigm (z) shown in the top panel 
of Figure 13 mini-halo clumping is strongly suppressed, and 



4.2 Suppression of H2 formation 

H2 formation and cooling in mini-halos will be suppressed 
relative to the no preheating case. This is because the fi- 
nite entropy of the gas allows it to resist compression; thus, 
the gas is at considerably lower density. CoUisional processes 
with cooling times tcooi oc 1/n will be suppressed relative to 
radiative processes which photo-dissociate H2. It has been 
recognized that gas profiles should have a central core rather 
than a cusp (e.g., Shapiro, Iliev & Raga (1999)) and that 
lower central densities will affect H2 chemistry (Haiman, 
Rees & Loeb 1996; Tegmark et al 1997). However, the fact 
that preheating can produce a much larger, lower density 
core than hitherto considered, and the subsequent implica- 
tions for H2 formation and cooling, has not been explored. 

It has been shown that H2 formation is enhanced in fos- 
sil HII regions due to increased electron fraction there (Ri- 
cotti, Gnedin & ShuU 2002) . Although H2 formation is very 
slow in the low density IGM, a high initial electron abun- 
dance will persist when weakly ionized gas is accreted onto 
non-linear structures, and aid H2 formation at that point 
(Oh 2000). We can conservatively take into account the H2 
enhancement by assuming the maximum initial abundance 
of H2, XH2 ~ 10~^; this hard upper limit is independent 
of density or ionization fraction and is due to 'freeze-out' 
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Figure 9. The evolution of the collapsed gas fraction in mini- 
halos (lower panel) and the gas clumping factor (upper panel) 
with redshift. The results are shown for no preheating (with 
^IGM = ^min)i well as constant values of entropy ii'iGM = 
1, lOeVcm'^ and the redshift-dependent entropy Kiqm(z) shown 
in the top panel of Figure 2. Both the global gas fraction in mini— 
halos and gas clumping are strongly suppressed for realistic values 
of the entropy floor. In particular, C ^ 1 for realistic values of 
KiQM, greatly reducing the photon budget required for rcioniza- 
tion. 



(Oh & Haiman (2002), hereafter OH02). The H2 abundance 
can only exceed this if three-body processes are important, 
which only takes place at very high density n > lO^cm"'^. 

In OH02, we showed that the minimum temperature 
Tmin a parcel of gas can cool to depends almost exclusively 
on tcooi/^diss oc Juv/n (see Figure 6 of OH02). This scal- 
ing behavior is only broken when the gas approaches high 
densities n > 10'*cm~^ (at this point, the cooling time be- 
comes independent of density and tcooi/idiss oc Juv)- How- 
ever, the latter regime is never reached in pre-heated gas, 
which is at much lower densities (n ~ 10*cm~^ corresponds 
to 5 ~ 6 X 10'' at z = 19) . In subsequent discussion, we shall 
assume the dependence of Tmin on Jw/n shown in Figure 
6 of OH02. 

Thus, in the presence of a radiation field Juv, the gas 
must be at a minimum density ricrit to cool down to a tem- 
perature Tmin- We have already calculated the maximum 
central density of gas in a halo given an entropy param- 
eter K, as in Figure |7| If the central density is less than 
the critical density, Uc < Ticrit, then none of the gas in the 
halo can cool down to Tmin- For a given entropy param- 
eter K, there is therefore a minimum radiation field Juv 
above which no gas can cool down to Tmin- We plot this in 
the top panel of Figure 1101 This plot is valid for all halos 
at all redshifts provided K and Juv are both rescaled ap- 
propriately (note that since n cx (1 -I- z)^, to keep Jw/n 
constant, Juv oc (1 -|- z)^). To get a sense of typical val- 



ues of Juv, the radiation field corresponding to ionizing 
photons per baryon in the universe is Juv ~ ^^71^.71;, (1 + 

z)VlO"^^ergs-^cm-^Hz-^sr-^ « lOn^ (tt)^ (where m 
is the comoving baryon number density, and ftp is the Planck 
constant). 

We see that an entropy floor greatly reduces the radi- 
ation field required to prevent gas cooling. For reasonable 
values of K, the reduction can be as much as four orders 
of magnitude. Thus, even if there is only a very weak radi- 
ation field, in the presence of an entropy floor effective H2 
cooling and star formation will be quenched. Note also that 
the cooling time exceeds the Hubble time for large values 
of K. This regime is reached only when gas densities are so 
low that even a very weak radiation field can dissociate H2: 
tcooi ~ tdiss ~ 2 X 10**(Juv/10-'*)-Vr. 

Without an entropy fioor, the density profiles of mini- 
halos at a given redshift are self-similar; thus, roughly com- 
parable fractions of their gas can cool down to Tmin and will 
be available as fuel for star formation. In the presence of an 
entropy floor, the fraction of gas with n > ricrit is greatly 
reduced, and much less gas can cool. Moreover, the self- 
similarity is broken: shallower potential wells (which have 
higher K) are more strongly affected by an entropy floor, and 
their central gas densities are much lower. Using our halo 
density profiles, we can calculate the fraction of gas above 
the critical density, fgas = M(,(> ncrit)/[(n(,/r2m)Mhaio]. We 
plot this in the lower panel of Fig llUI At some critical value 
of K, the fraction of gas which can cool plummets dramati- 
cally. This value of K corresponds roughly to when the low 
density core falls below the critical density. 

Note that we assume static rather than evolving density 
profiles. In reality, as gas cools, densities should increase as 
the gas loses pressure support, making the gas more resistant 
to H2 photo-dissociation (Machacek et al. 2001). However, 
the dynamical time tdyn ~ 10^(n/l cm~^)~^''^yr is typically 
much longer than the photo-dissociation time tdiss ~ 2 x 
10* Juy/J^Jj,[jyr (where the factor /siiidd takes into account 
self-shielding; see discussion in section 4.2.3 of OH02, and 
references therein). So our use of the initial density profile 
to estimate timescales should be approximately valid. 

We can also compute the critical radiation field re- 
quired to suppress H2 cooling down to some temperature 
T as a function of redshift. We plot this in Fig. 1111 for a 
Tvir = 8000K halo; such a halo represents an upper mass 
limit to the mini-halo population. Smaller halos, which have 
lower central densities for a given entropy floor, will require 
a weaker radiation field for H2 cooling to be suppressed. This 
is shown for four different entropy floors: Kigm{z) as shown 
in Fig. 2, no preheating A"igm ~ Jfmin, and fixed entropy 
floors KjGM = 1, lOeVcm^. Gas cooling in fossil H2 regions 
is very easily suppressed; the radiation field required is 2-4 
orders of magnitude lower than required in the no preheating 
case. 

One possible caveat is that since fossils are likely to be 
metal-enriched, gas cooling is dominated by metal line cool- 
ing, rather than H2 cooling. We argue that this is unlikely, 
for two reasons: 1) Metal line cooling only dominates H2 
cooling above a critical metallicity Z > 10~^ Zq (Bromm et 
al 2001; Hellsten & Lin 1997). This level of metal enrich- 
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Figure 10. Top panel: The critical radiation field Juv required 
to suppress gas cooling below some temperature T, as a function 
of entropy parameter K. An entropy floor greatly reduces the 
radiation field required to photo-dissociate H2, by several orders 
of magnitude. Note that Juv oc (1 + 2)^; the curves shown are 
for a halo at z = 14. Bottom panel: The fraction of gas in a 
halo able to cool below temperature T, for Jjjv = 10~^ (lines), 
and Juv = lO"'^ (dots connected by lines), at z = 14. Beyond a 
critical value of K, the cooled gas fraction plummets dramatically. 

ment corresponds to ~ 1/^^ ionizing photon per baryon in 
the universe (where fz is the volume filhng factor of met- 
als). Thus, metal line cooling becomes significant only at late 
times. 2) After an ionizing source turns off and explodes 
as a supernova, the metal-polluted region is much smaller 
than the fossil HII region (Madau, Ferrara & Rees 2001). 
Most of the volume is therefore still of pristine composi- 
tion, and undergoes the entropy floor suppression we have 
described. The metal-polluted region lies close to the high 
density peak where the very first stars formed, where in any 
case, more massive halos Tvir > iCK will collapse. To sum- 
marize: metal line cooling is therefore unlikely to spoil our 
assumption of adiabaticity in most of the volume of the fos- 
sil HII region. The exception is at the most highly biased 
density peaks, where metals can be thought of as effectively 
increasing the star formation efficiency. 

The net result is that entropy injection greatly boosts 
the negative feedback from early star formation: the entropy 
floor in reionized regions results in low density cores in the 
center of halos, in which H2 is easily photo-dissociated by a 
weak external UV radiation field. 

4.3 Negative Feedback from X-rays 

It is often argued that X-rays boost H2 production and cool- 
ing in mini-halos, by penetrating deep into the dense core 
and increasing the free electron fraction, which is critical for 
gas phase H2 production. Thus, X-rays are thought to exert 



Figure 11. The critical radiation field Juv required to suppress 
gas cooling below T = lOOOK, as a function of redshift z. The 
curves are shown for a Tvir = 8000K halo, which is an upper 
limit in mass to the mini— halo population: smaller mini— halos 
will require a weaker radiation field for cooling to be suppressed. 
The curves axe shown for four different entropy floors: no pre- 
heating K'iGM = Jfmin, an evolving entropy Kiqm_{z) as shown 
in Fig. 2, and flxed entropy floors Kiqm = 1, lOeVcm^. Points 
are also shown for 2ct and 3cr halos at each redshift, assuming 
KiGuiz) (the points end when Tvir < lOOOK and Tvir > 8000K 
respectively). Typical entropy floors at each redshift reduce the 
critical radiation field required for H2 suppression by 2 — 4 orders 
of magnitude, compared to the no preheating case. 



a positive feedback effect, counter-acting photo-dissociation 
by UV radiation (Haiman, Rees & Loeb 1996; Haiman, Abel 
& Rees 2000). We show here by fairly general arguments 
that X-rays in fact exert a strong negative feedback effect, 
due to the entropy that they inject into the IGM. This pre- 
vents gas from compressing to sufficiently high density to 
produce H2 and cool efficiently, and far outweighs any pos- 
itive feedback from increasing the free electron fraction. 

For the sake of dofiniteness, we consider an intrinsic 
source spectrum of the form 

J = Juv (^^) ^ hiy < 13.6eV (22) 
Jx ^ 13.6eV <hiy< lOkeV 

where huL = 13.6eV and Jvv/Jx ~ /break = 1 — 100 is 
the spectral break at the hydrogen Lyman edge; note that 
fbreak > 1 always. A spectrum of the form Ji, oc is chax- 
acteristic of the mean spectra of quasars (Elvis et al 1994), 
as well as inverse Compton emission from high-redshift SN 
(Oh 2001). It is distinctive in that vJi,=consi, i.e., there 
is equal power per logarithmic interval (though our argu- 
ments can be generalized to other spectra). Note that the 
actual ionizing spectrum at any given point in the IGM will 
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be much harder, since photoelectric absorption hardens the 
spectrum away from the source. 

What is the heating associated with this radiation field? 
As we argued in H2.2l above. all energy in the X-ray radiation 
field with El < E < Ethick (as defined in fj^l will be absorbed 
by the IGM, since the universe is optically thick at these 
frequencies. The radiation field may be subdivided into two 
components. The "mean field" consists of photons with E > 
-Eovoriap, where ^overlap is defined by: 

Amfp(£) > 7^som?c for E > -Bovcrlap (23) 

i.e., photons with E > iSovcriap (typically, iJovoriap ~ lOOeV) 
have a mean free path greater than the mean separation be- 
tween sources, so that a homogeneous ionizing background 
is established. The "fiuctuating field" consists of photons 
with E < -Eoveriap; this component is dominated by radi- 
ation from a single source, and is subject to large Poisson 
fluctuations. We are interested in the heating due to the 
'mean field', which can lead to an entropy fioor even outside 
the HIT regions of ionizing sources. It is: 

-Eheat = /heat / due^, ^ /heat 7^ Jx In f (24) 

where /heat is the fraction of energy of the hot photo-electron 
created by an X-ray which goes into heating. We have con- 
servatively set ~ {4:nJ^)/lH', this underestimates for a 
given J^. For the spectrum we have chosen, the result is only 
logarithmically sensitive to the integration bounds; we set 
In (fthick/^'thin) ~ 2.3. For steeper spectra, oc z/~" where 
Q > 1, the lower integration limit becomes important. 

Over a Hubble time tH, the internal energy density U — 
nksT ~ -Bhcatiii of the gas becomes: 

C/«/heat — Jx/^Llnf^^) (25) 

This make sense: since some fraction /heat of the radiation 
field goes directly into heat, /hcat47r(i/Jiy) ~ C/c. Thus, X- 
rays will heat the IGM to a temperature: 

There is therefore a direct relation between the X-ray radi- 
ation field and the entropy fioor: 

A-,™.l..eVc„=(^)(^)(i±f)-, 

As we argued in H2. 21 above, the IGM evolves adiabati- 
cally for Xe < 0.1[(1 -I- z)/15f-^. What level of the radiation 
field does this correspond to? For Xe < 0.1, about /ion ~ 1/3 
of the injected energy goes toward performing ionizations; 
the number of photo-ionizations is simply -Ephoton/(37eV). 
Thus, we obtain: 

..^2.3xl0-(^)(^) (28) 

and so for Jx < 5, the IGM evolves adiabatically (note 
also that since Uec ~ 10te(xe/10~^)"^(T/100K)-'' '^([l -f 
z]/15)~^'^ , we can ignore recombinations). 

Given the entropy fioor from equation 1271 we can now 



calculate density profiles of halos. In particular, we can cal- 
culate the maximum central density of a halo from the en- 
tropy floor induced by a given X-ray background. As be- 
fore, from the central density of a halo we can calculate the 
critical level of the LW flux Juv required to suppress H2 
cooling down to some temperature Tmin. This critical flux is 
shown as a function of /x = Jx/ Jvv in Fig |12l for a halo of 
Tvir = 5000K at 2; = 19. The curves may be simply rescaled 
for halos of different mass at different redshifts by trans- 
lating the X-axis: fx oc Tvir(l + z)^. The effect of X-rays 
increases for less massive halos (which are less able to with- 
stand an entropy floor) at lower redshifts (when gas densities 
are lower) . Also shown are the maximum value of the X-ray 
background /x ~ 1 for the hardest sources, as well as the 
typical minimum value of the X-ray background required 
for positive feedback to be effective /x ~ 0.1 (Haiman, Abel 
& Rees 2000). The shaded region therefore denotes the X- 
ray fluxes thought to produce positive feedback effects. We 
see for these large values of the X-ray background. X-rays 
cannot exert a positive feedback effect. The required X-ray 
background is so large that it would produce an unaccept- 
ably large entropy floor, which prevents the formation of 
dense regions where H2 can withstand photodissociation. In- 
stead, X-rays produce an over-all negative feedback effect. 
In fact, values of /jc « 0.1 — 1 reduce by almost two orders 
of magnitude the value of the UV background required to 
photo-dissociate H2 and prevent it from cooling down to low 
temperatures. 

The only way of circumventing these difficulties is if 
X-rays turn on suddenly at some epoch. Halos which have 
already collapsed will be unaffected by the entropy floor, 
but will undergo the usual positive feedback from X-rays. 
However, the subsequent generation of mini-halos will be 
suppressed by the entropy floor. 



4.4 21cm Observational Signatures 

At present, the only observational probes proposed for small- 
scale structure at high redshift such as mini-halos are 21cm 
observations with future radio telescopes such as the Square 
Kilometer Array (SKA)^ and the Low Frequency Array (LO- 
FAR)"^. It has been proposed that mini-halos could be de- 
tected statistically in emission (Iliev et al 2002a; Iliev et al 
2002b) and individually in absorption along the line of sight 
to a high-redshift radio source (Furlanetto & Loeb 2002). 
If an entropy floor exists, such signatures will be strongly 
suppressed. 21cm signatures of mini-halos therefore provide 
an interesting indirect probe of reionization history: if mini- 
halos are detected in large numbers at a given redshift, it 
would imply that the flUing factor of fossil HII regions is 
still small at that epoch. In particular, a comoving patch 
where mini-halos are seen has likely never been ionized be- 
fore (which is a much stronger constraint than merely being 
neutral at the observed epoch). This could prove a very use- 
ful probe of the topology of reionization. 

^ see http:/ /www. nfra.nl/skai" 
^ see http: / /www.astron.nl/lofar| 
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Figure 12. The UV background Juv required to suppress cooling 
down to some temperature Tmi^, as a function of /x = Jx/J\jv- 
The curves plotted are for a Tvir = 5000K halo at z = 19; the 
results may be rescaled to other halos by rescaling the x-axis: 
fx T^iii^ + z)^- Also shown as dashed and dotted lines are 
the maximum value of fx ^ 1 in the hardest sources, and the 
minimum typical value of /x ~ 0.1 required for X-rays to ex- 
ert a positive feedback effect; the shaded region therefore shows 
the possible region whereby X-rays were assumed exert a pos- 
itive feedback effect. Instead, we see that such a strong X-ray 
background exerts a negative feedback effect. By inducing a low 
density core in the halo, an X-ray background in the required 
zone reduces by 1 — 2 orders of magnitude the minimum UV flux 
required to suppress cooling, compared to the case where fx —* 0. 



Mini-halos are too faint to be seen individually in emis- 
sion, and can only be detected statistically through bright- 
ness temperature fluctuations. Provided Ts 3> Tcmb (where 
Ts is the spin temperature), the 21cm flux S is independent 
of the spin temperature, and depends only on the HI mass, 
S oc Mhi- Thus, S'(halos)/S'(IGM) is simply equal to the 
collapsed gas mass fraction in mini-halos, which is always 
less than unity. The mini-halo signal is much smaller if an 
entropy floor exists, since the collapsed gas fraction declines 
rapidly: the bottom panel of Figure|3|can simply be read off 
as S(halos)/S(IGM). Thus, the IGM dominates 21cm emis- 
sion. Mini-halos dominate only when: (i) the IGM spin tem- 
perature has not yet decoupled from the CMB (the critical 
thermalization radiation fleld J21 such that Ts ^ Ta ~ Tk, 
where Ta is the color temperature of the radiation field, is 
J21 « 5([l-|-«]/10) (Madau, Meiskin & Rees 1997)) (n) the 
filling factor of fossil HII regions is small. The high WMAP 
optical depth, which implies significant reionization at high 
redshift, pushes the latter constraint to high redshift and 
thus low radio frequencies Uobs ~ 93([1 + 2]/15)MHz, mak- 
ing this a very challenging observation. 

We now turn to absorption signatures. How does the 



mini-halo 21cm absorption signature scale with the IGM 
entropy? The IGM 21cm optical depth is: 



/ Tcmb{z) 



10 



1/2 



xni. 



(29) 



See Madau, Meiskin & Rees (1997) for expressions for the 
spin temperature Ts{J2i,nHi,TK); we do not reproduce 
them here. We use a fit to the results of Allison & Dal- 
garno (1969) for the coUisional coupling coefficient Cio for 
Tk < lOOOK, with a Cio oc T~°-^^ extrapolation for higher 
temperatures. The mini-halo 21cm optical depth along an 
impact parameter a to the halo center is (Furlanetto & Loeb 
2002): 



Thalo(a, V) 



Shpc^Aw 



dR 



nHi{r) 
Ts{r)^b{r) 



exp 



viyf 
&2(r) 



(30) 



where Aio — 2.85 x lO^^^s^^ is the spontaneous emission co- 
efficient, h{r)'^ = IksTK /rup is the gas Doppler parameter, 
= + a^, and v{i') = c{y — Vo)/vo)- Unlike Furlanetto 
& Loeb (2002), we do not attempt to model the velocity 
field of the infall region, which we ignore; the total observed 
optical depth is taken to be rtot = Thaio + tigm- The veloc- 
ity field in the preheated gas will differ significantly from 
the Bertschinger self-similar solution they assume, since gas 
pressure retards accretion; the infall region optical depth 
will be smaller because of the reduced gas column. Since 
even in their case the contribution of the infall region to 
21cm absorption is small (see their Fig. 2, where the equiv- 
alent width drops drastically for ct > rshock), our neglect is 
justified. 

In Fig. 1131 we show the optical depth as a function of 
observed frequency, for different levels of the entropy param- 
eter K. The values of K shown correspond to IGM temper- 
atures of TiGM = 70, 700, 7000K; since fossil HII regions at 
z=10 cannot significantly cool below 7000K (see Fig.0, the 
value of = 50 is most appropriate. We see that the opti- 
cal depth contribution from mini-halos falls drastically with 
increasing K. Also shown is the optical depth for J21 = 10. 
The radiation field drives Ts Tk, reducing the 21cm op- 
tical depth r oc Tg^ . 

Since tigm also depends on K, the most relevant quan- 
tity for detecting a mini-halo is the observed equivalent 
width: 



(31) 



which measures the 21cm absorption due to the mini-halo in 
excess of that due to the IGM. It falls rapidly with K and is 
extremely small for the most likely value oi K = 50. Figures 
1131 and IT^ can be compared with Fig. 1 and 2 of Furlanetto 
& Loeb (2002) (although note that in Fig. 1131 we plot r as 
a function of observed rather than intrinsic frequency) . An 
entropy fioor greatly reduces the cross-sectional area over 
which an observable absorption signal may be detected. If 
one detects significant mini-halo absorption over a large co- 
moving patch along the line of sight to a radio source, one 
can place an upper limit on the entropy fioor there. This in 
turn places an upper limit on the high-redshift X-ray back- 
ground (from equation IHJ , as well as a lower limit on the 
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Figure 13. Mini-halo optical depth profiles as a function of ob- 
served frequency for a 2 X 10^Mq/i~^ halo at z = 10, all for impact 
parameters of a = O.Srvir- The values of K = 0.5,5,50 shown 
correspond to IGM temperatures of Tjqm = 70, 700, 7000K; since 
fossil HII regions at z=10 cannot significantly cool below 7000K, 
the value of = 50 is most appropriate. Solid lines correspond 
to J21 = 0; dashed lines correspond to J21 = 10. The entropy 
floor makes the mini-halo contribution to the 21cm absorption 
undetectable. 



redshift at which that patch was first reionized (since fossils 
from higher redshift have lower entropy, as in Fig|5J. 

The chief driver of high-redshift 21cm proposals has 
of course been to observe the IGM itself in 21cm emission 
(Madau, Meiskin & Rees 1997; Tozzi et al 2000; Ciardi 
& Madau 2003). We note that if the filling factor of fossil 
HII regions is large, conditions are very favorable for such 
observations. Provided Ts —> Tk in the IGM, the only con- 
dition for the IGM to be seen in 21cm emission against the 
CMB is for Tk S> TcAfs(-Zobs)- This is certainly satisfied in 
fossil HII regions. It was previously thought that recoil heat- 
ing from Lya photons could heat neutral regions (Madau, 
Meiskin & Rees 1997), but a detailed calculation (Chen & 
Miralda-Escude 2003) shows the heating to be insufficient. 
In that case, a high-redshift X-ray background would be re- 
quired to heat neutral regions. Such concerns are moot if a 
period of early reionization took place (as seems to be indi- 
cated by WMAP observations), followed by recombination: 
large tracts of warm, largely neutral gas would exist. Note, 
however, that brightness temperature fluctuations can only 
be detected if the contribution from unresolved radio point 
sources (Tozzi et al 2000; Oh & Mack 2003) can be suc- 
cessfully removed by using spectral structure in frequency 
space. 



Figure 14. The observed 21cm absorption equivalent width of a 
2 X IO'^Mq/i-i halo at z = 10, for J21 = 0, 1, 10, as a function of 
K. Solid lines are for an impact parameter of a = O.Srvir, while 
dashed lines are for a = 0.7r^\^. The equivalent width falls rapidly 
as the entropy floor rises; for the expected value of if = 50 at this 
redshift, equivalent widths are very small. 

5 EFFECTS ON GLOBAL REIONIZATION 
SCENARIOS 

5.1 General Considerations 

We next address the effect of the suppression of mini-halo 
formation on the global reionization history. The importance 
of the effect depends on (1) the synchronization of the for- 
mation of mini-halos tsync relative to the typical lifetime 
tion of the ionizing source, which determines the fraction of 
mini-halos subjected to feedback, as well as (2) on the re- 
combination time ircc relative to tsync, which determines how 
long the ionized regions last and thus the fraction of active 
(as opposed to fossil) ionized regions at any given time. 

In the limit of taync <C iion, the fossil HII regions would 
appear only after all the mini-halos had already formed, 
and hence the feedback would have no effect on the total 
amount of reionization by mini-halos. In the absence of any 
other feedback effects, the IGM could then, in principle, be 
fully reionized by mini-halos, given a high enough ionizing 
photon production efficiency (although would subsequently 
recombine) . 

In the opposite limit, tsync S> tion, the contribution of 
mini-halos to reionization will be strongly suppressed. We 
can obtain a rough estimate for the maximum fraction of 
the IGM that can be ionized by the mini-halos. During the 
lifetime of its resident ionizing source, each mini-halo pro- 
duces an ionized volume Vhii, which will correspond to the 
total number N-y of ionizing photons injected into the IGM, 
i.e. nV ~ A'^^, where n is the mean hydrogen density. To a 
good approximation, recombinations can be ignored in this 
phase, since massive stars have lifetimes shorter than the 
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recombination time, Uec « 3 x 10^[C(1 + 2)/25]~^ years (see 
Fig. llSl below'). On the other hand, the recombination time is 
shorter than the Hubble time, froc/fhub ~ [C(l + 
so in general, at high redshift z > 10, and for clumping 
factors C > 1, the ionized volume will recombine once its 
driving source turns off. 

As argued above, for the purpose of the mini-halo sup- 
pression, we may nevertheless simply imagine that the HII 
region never recombines. Ignoring recombinations inside the 
ionized region, this is equivalent to setting the size of the 
fossil HII region to be the maximum (comoving) size of the 
active HII region, reached at the time when the ionizing 
source turns off. 

The evolutions of the radii of the active and fossil HII re- 
gions are illustrated in Figure lT^ for a single 10^ Mq metal- 
free star that turns on at z = 25. We follow the expansion of 
the ionization front _Rs by solving the standard differential 
equation, taking into account ionizations, recombinations, 
and the Hubble expansion (e.g. Shapiro & Giroux 1987; Cen 
& Haiman 2000). We assume a constant clumping factor of 
C = 1 (upper solid curve) or C = 10 (lower solid curve). We 
assume further that the metal free star emits 40, 000 ionizing 
photons per stellar proton at a constant rate for ~ 3 x 10® 
years (Tumlinson & ShuU 2000; Bromm, Kudritzki & Loeb 
2001; Schaerer 2002). Once the ionizing source turns off, the 
solid curves show the formal solution for the evolution of the 
ionization front. * The dashed curves show the size of the 
fossil HII region. 

We can next consider the evolution of the volume filling 
factor of the fossil HII regions from the ensemble of mini- 
halos, and define the epoch Z{ when the filling factor of these 
fossils reaches unity. As argued above, no new mini-halos 
can form at z < Zf . At this epoch, the global ionized fraction 
will be smaller than unity, because each fossil HII region has 
already partially recombined. Thus the global ionized frac- 
tion will be ~ exp(— taync/iroc), assuming that the sources 
typically turned off a time isync ago, and free is the aver- 
age recombination time over the interval tsync. For example, 
assuming that the fossil HII regions overlap at Zf — 20, 
with the typical sources born at z = 25, and the recombi- 
nation time evaluated at z = 22.5 (with C — 1), we find 
tsync ~ tree ~ 5 X 10^ ycars, and a maximal ionized fraction 
of ~ 30%. 

5.2 Models for the Reionization History 

To refine the above considerations, we next compute the evo- 
lution of the global ionized fraction using a semi-analytical 
model adopted from Haiman & Holder (2003, hereafter 
HH03). For technical details, the reader is referred to that 
paper. Here we only briefiy summarize the main ingredients 
of the model, and describe the modifications we made to 

* This formal solution predicts a shrinking of the ionization front. 
In reality, once the source turns off, recombinations throughout 
the ionized region will decrease the ionized fraction uniformly 
by the factor [_Rs/_Rg maxl'^i rather than reduce the size of a 
highly ionized volume. However, for our purposes of computing 
the global ionized fraction, the two interpretations are equivalent. 
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Figure 15. The evolution of comoving radius of the ionized region 
around a single 10^ Mq metal-free star that turns on at 2: = 
25, and is assumed to emit 40, 000 ionizing photons in ~ 3 X 
10^ years. The upper and lower solid curves show solutions that 
include recombinations with C = 1 and C = 10. The ionized 
volume shrinks (or equivalently, the ionized fraction is decreased 
by the factor [-Rs/Rs.max]"^) after the source turns off at 2 ~ 24. 
The corresponding dashed curves show the radii of the fossil HII 
regions; these are assumed to stall at the maximum radius of the 
HII sphere and stay constant thereafter. 



include the suppression of the formation of new mini-halos 
inside fossil HII regions. 

In this model, we follow the volume filling fraction Fmi 
of ionized regions, assuming that discrete ionized Stromgren 
spheres are being driven into the IGM by ionizing sources lo- 
cated in dark matter halos. For simplicity, we ignore Helium 
in this work. The dark matter halo mass function is adopted 
from Jenkins et al. (2001). We consider two distinct types 
of ionizing sources, located in halos with virial temperatures 
of 10^ K ^ Tvir <, W^K (mini-halos) or Tvir ^ 10*K (large 
halos). We assume that in mini-halos, a fraction /, — 0.005 
of the baryons turn into metal free stars (Abel, Bryan & 
Norman 2000, 2002; Bromm, Coppi & Larson 1999, 2002), 
which produce A'^^ = 40, 000 ionizing photons per baryon, 
and fcsc ~ 100% of these photons escape into the IGM. The 
product e = /*/escA.^ — 200 determines the size of the ion- 
ized regions around mini-halos. For large halos, we adopt 
/, = 0.1, /esc = 20%, and = 4,000, for a total efficiency 
e — 80. In contrast with HH03, we here do not distinguish 
halos with 10* K < Tvi^ < 2 x 10^ K from Tvir > 2 x 10^ K 
(see Dijkstra et al. 2003 for the reason why this distinction 
is likely unimportant). 

As in HH03, we then follow the evolution of the global 
ionized fraction Fhii(2) by summing the ionized volumes 
surrounding the individual dark matter halos. We exclude 
the formation of mini-halos in active ionized regions (but 
assume larger halos are impervious to photoionization feed- 
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back). The evolution of Fhii in this model is shown by the 
upper thin solid curve in Figure^] (with C = 1 in the upper 
panel, and C = 10 in the lower panel). As the figure shows, 
in the absence of other feedback effects (see discussion in 
Haiman, Abel & Rees 2000 and Haiman 2003 for other feed- 
back effects), the mini-halos could ionize the IGM in full by 
2 ~ 14 (C = 1), or nearly fully by 2 ~ 10 (if clumping is 
assumed to be more significant, C — 10). 

We next compute _Fhii in the same model, except that 
we exclude the formation of new mini-halos in fossil HH re- 
gions. This is easily accomplished in practice: the suppres- 
sion factor (1 — Fhii) for the formation rate of mini-halos is 
replaced by a factor (1 — -F/m). Here the fossil filling factor 
Fhjj is computed the same way as the ionized fraction -Fhii, 
except that the individual ionized regions are assumed to fol- 
low the dashed curves from Figure \TE\ rather than the solid 
curves. In Figure ll(il we show the evolution of the volume 
filling factor of the fossil HH regions (dashed curves) , 
as well as the reionization history Fhii{z) as the thick solid 
curves. Finally, for reference, we show Fhii(«) with mini- 
halos completely excluded (e — 0) as the lower thin solid 
curves. 

As apparent from Figure [T^ the exclusion of new mini- 
halo formation from the fossil HH regions causes a significant 
suppression of the total ionized fraction that can be reached 
by mini-halos. Under the rather optimistic set of assump- 
tions described by the thick solid curve in the upper panel, 
the maximum ionized fraction that can be reached is ~ 40%. 
In reality, clumping is unlikely to be unity in the immediate 
vicinity ( <, 100 kpc) of the ionizing sources (Haiman, Abel 
& Madau 2000). 

The total electron scattering optical depth attributable 
to mini-halos (the appropriately weighted integral between 
the thick solid curve and the lower thin solid curve) is 
T = 0.07 and r = 0.014 in the C = 1 and C = 10 cases, re- 
spectively. This makes it unlikely that mini-halos can fully 
account for the large optical depth t — 0.17 measured by 
WMAP. Note the thick curve in the upper panel of Figure lTSl 
has a total r — 0.2 (r = 0.07 attributable to minihalos, and 
T = 0.13 to larger halos), so that it is consistent with the 
WMAP measurement. Note that raising the efficiencies in 
mini-halos would not increase the optical depth attributed 
to minihalos, since feedback then would set in earlier (we 
have explicitly verified that r is approximately independent 
of efficiencies over a range of multiplicative factors 0.1-10 
for e) . Finally, note that the suppression considered here and 
shown in Figure^] provides a negative feedback in addition 
to the negative feedback expected from H2 photodissociation 
(Haiman, Rees & Loeb 1997; Haiman, Abel & Rees 2000). 
In fact, as discussed above, preheating amplifies the effect 
of the H2-photodissociative negative feedback. 

Because an entropy floor essentially eliminates gas 
clumping due to mini-halos (as shown in section 14. H , the 
reionization history is likely to more closely approximate 
the C=l case than the C=10 case until low redshifts z < 10 
(when gas clumping due to larger structures predominates). 
An entropy floor therefore has two countervailing effects on 
reionization: by suppressing star formation in mini-halos, it 
reduces the comoving emissivity. However, by reducing gas 




z 



Figure 16. The evolution of the volume filling fraction -Fhii of 
ionized hydrogen in reionization models with different assump- 
tions about feedback on mini— halo formation. The upper and 
lower panels assume constant clumping factors of C = 1 and 
C = 10, respectively. In both panels, the dashed curves show the 
volume filling fraction of the fossil HIT regions. The thick solid 
curves show the evolution of -Fhii assuming that mini-halos can- 
not form inside fossil HIT regions. The upper solid curves assume 
that mini-halos are only excluded from forming inside active ion- 
ized regions. Finally, the lower solid curves ignore the contribution 
to Fhii from any mini-halos. 

clumping, it also reduces the photon budget required for 
reionization. It is interesting to note that that the evolution 
of the filling factor is non-monotonic for the C=l case: feed- 
back due to early reionization naturally produces a bump 
in the comoving emissivity, and a pause ensues before larger 
halos (which can resist feedback) collapse. This is similar to 
the reionization histories derived by (Cen 2003; Wyithe & 
Loeb 2003), but is regulated by feedback from the entropy 
fioor rather than a Pop HI to Pop II transition due to a 
universal metallicity increase. 



6 CONCLUSIONS 

In this paper, motivated by the WM-AP results, we have con- 
sidered the feedback effect of early reionization/preheating 
on structure formation. This feedback effect is inevitable in 
any reionization scenario in which star formation throughout 
the universe is not completely synchronized. Our principal 
conclusions are as follows: 

1. Fossil HH regions have a residual entropy floor after 
recombination and Compton cooling which is higher than 
the shock entropy for mini-halos (Tvir < lO^K); thus, such 
halos accrete gas isentropically. The IGM entropy depends 
primarily on the redshift and only weakly on the overden- 
sity (5; it is thus largely independent of the details of struc- 
ture formation. For this reason, and also because it is con- 
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served during adiabatic accretion or Hubble expansion, the 
gas entropy is a more fundamental variable to track than the 
temperature. We provide a simple analytic formula for the 
temperature (and hence entropy K = T/n^^^), in equations 
1^ and 1211 An early X-ray background would also heat the 
entire IGM to similarly high adiabats. The entropy floor due 
to the latter would be much more spatially uniform. 

2. We apply the entropy formalism used to calculate 
the effect of preheating on low-redshift galaxy clusters to 
the high-redshift minihalos. We obtain detailed gas density 
and pressure profiles, which we use to calculate the impact 
of preheating on the central density, accreted gas fraction, 
gas clumping factor, and mini-halo baryonic mass function. 

3. These quantities can then be used to calculate global 
effects of preheating. The collapsed gas fraction in minihalos 
falls by ~ 2 orders of magnitude, while the gas clumping 
factor falls to C ^ 1, as for a uniform IGM. An entropy floor 
reduces the photon budget required for full reionization by 
about an order of magnitude, by reducing gas clumping and 
eliminating the need for mini-halos to be photo-evaporated 
before reionization can be completed. 

4. However, an entropy floor does not necessarily pro- 
mote early reionization: it also sharply reduces the comoving 
emissivity. By reducing the central gas densities in mini- 
halos, preheating impedes H2 formation and cooling, and 
reduces the critical UV background required for H2 sup- 
pression by 2-4 orders of magnitude. Thus, once a comoving 
patch of the IGM is reionized, no subsequent star formation 
in mini-halos can take place in that volume. The patch can 
only be reionized by more massive halos r > 10*K, which 
can undergo atomic cooling. By furnishing an entropy floor, 
X-rays also suppress H2 formation. Thus, contrary to con- 
ventional wisdom. X-rays provide negative rather than pos- 
itive feedback for early star formation. 

5. Mini-halos will not be observable in 21cm emis- 
sion/absorption in fossil HII regions. Thus, 21cm observa- 
tions provide an unusual probe of the topology of reioniza- 
tion: mini-halos trace out regions of the IGM which have 
never been ionized. If mini-halos are seen in large numbers, 
this places an upper limit on the filling factor of fossil HII 
regions and the X-ray background at that redshift. 

6. We have computed the reionization histories as in 
HH03, but taking the feedback effect of an entropy floor 
(and reduction of gas clumping) into account. The strong 
feedback in fossil HII regions imply that HII fronts at high 
redshift never overlap, and global reionization at high red- 
shift does not occur. This limits the contribution of mini- 
halos to the reionization optical depth r < 0.07, almost in- 
dependent of star formation efficiency in mini-halos (if star 
formation is more efficient, feedback sets in earlier). Thus, 
the bulk of the optical depth observed by WMAP must come 
from more massive objects. 

7. Strikingly, we obtain a double-peaked reionization 
history: an early peak in which the universe is filled with fos- 
sil HII regions, followed by a pause before more massive ha- 
los collapse which finally fully reionize the universe. This is 
similar to 'double reionization' scenarios computed by other 
authors (e.g. Cen (2003); Wyithe & Loeb (2003)), but one 
in which the the comoving emissivity is regulated by gas en- 



tropy, rather than a Pop HI to Pop II transition due to a 
universal metallicity increase. This last point deserves ad- 
ditional comment. A metallicity- regulated evolution of the 
emissivity requires that metal pollution is fairly spatially 
uniform. However, different parts of the IGM likely undergo 
the Pop HI to Pop II transition at different epochs. Further- 
more, an increase in metallicity does not necessarily result 
in a drop in the overall emissivity: metal line cooling likely 
results in a larger star formation efficiency in halos, since 
metals are not subject to internal UV photo-dissociation, 
unlike H2. The factor of ~10 drop in the HI ionizing emis- 
sivity per stellar baryon could be outweighed by the increase 
in the total mass of stars formed. In comparison, we argue 
that an entropy-regulated transition is inevitable, and there- 
fore more robust. This is an important conclusion, given the 
possibility that future CMB polarization studies will be able 
to distinguish among different reionization histories (HH03; 
Holder et al 2003). 

In this semi-analytic study, we essentially assumed a 
single value of the entropy floor at each redshift, for all ha- 
los. In reality, as we discussed above, there should be large 
fluctuations in entropy, depending on the topology /history 
of reionization and the accretion/merger history of halos. It 
would be interesting to study these effects in detail in three 
dimensional numerical simulations. To begin with, it would 
be interesting to re-run existing simulations of the collapse 
of gas in mini-halos, but adding an entropy floor in the ini- 
tial conditions. Such studies could quantify more precisely 
the effects of an entropy floor in suppressing H2 formation 
and cooling in mini-halos. In fact, since densities are lower 
and gas cooling is reduced, it should be a computationally 
more tractable problem. With the aid of larger volumes, 
several global issues mentioned above could be addressed: 
for instance, a global, self-consistent study of feedback as in 
Machacek, Bryan & Abel (2001), but including the effects 
of both the entropy floor and UV feedback. Because of its 
strong effect on the evolution of the comoving emissivity, gas 
entropy acts as a self-regulating mechanism which likely has 
a strong influence in controlling the progress of reionization. 
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APPENDIX I: ANALYTIC EXPRESSION FOR 
THE FINAL TEMPERATURE 

It is useful to have an approximate analytic expression for 
the final temperature a parcel of gas cools down to, given the 
initial temperature Ti, redshift z, overdensity 5 and length 
of time spent cooling t. This allows one to quickly estimate 
the effect of early reionization in different situations with- 
out evolving the full chemistry code. We develop such an 
expression in this Appendix. 

At the redshifts and overdensities of interest, Comp- 
ton cooling dominates by far. The most relevant physics is 
Compton cooling and hydrogen recombination: 



T 



8 aT^ar 



3 rrieC 1 + X, 



■{T-T^) 



(32) 



Note that (T - T^) « T since T > T^, and a oc T"" '^ in 
the temperature range of interest. Assuming the gas is fully 
ionized = 1 at some initial temperature Ti, we obtain the 
analytic solution: 
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(33) 



(34) 



This gives temperature as a function of ionization fraction 
Xe- We therefore need to know the final ionization fraction. 
How can we estimate it? If the gas recombines isothermally 
at temperature T' , the ionization fraction is given by: 



Xe{t) = 



1 + {t/treciT')) 



(35) 



where Xo = 1 is the initial ionization fraction. If the gas cools 
as it recombines, substituting the instantaneous tempera- 
ture T' into the expression, would overestimate the speed 
of recombination and underestimate x^. We find that if we 
substitute T' = 2T and substitute this into equation 1331 
the solution of this non-linear equation for T/ is remarkably 
close to the full non-equilibrium solution (see points in Fig- 
ure 0. We have also verified that we obtain fairly accurate 
results for Xe{t). Equations I33l and l35l thus give Tf{z, 5, Ti,t). 
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Fossil HII Regions: Self-Limiting Star Formation at High Redshift 

Our neglect of recombination line cooling fails in high 
density regions. This leads to at most a factor ~ 2 error in 
the final temperature, since recombination line cooling can 
cool gas down to at most ~ 5000K. The entire discussion 
assumes that H2 formation and cooling is not competitive 
with Compton cooling, which is generally true in low density 
regions. 
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